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Abstract 


This study examines the importance of horizontal photon transport effects, 
which are not considered in the 1-D calculations of solar radiative heating used by 
most atmospheric dynamical models. In particular, the paper analyzes the difference 
between 2-D and 1-D radiative calculations for 2-D vertical cross-sections of clouds 
that were observed at three sites over 2- to 3-year periods. The results show that 
2-D effects increase multiyear 24-hour average total solar absorption by about 
4.1 W/m 2 , 1.2 W/m 2 ' and 0.3 W/m 2 at a tropical, mid-latitude, and arctic site, 
respectively. However, 2-D effects are often much larger than these average values, 
especially for high sun and for convective clouds. The results also reveal a somewhat 
unexpected behavior, that horizontal photon transport often enhances solar heating 
even for oblique sun. These findings underscore the need for fast radiation 
calculation methods that can allow atmospheric dynamical simulations to consider 
the inherently multidimensional nature of shortwave radiative processes. 
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1. Introduction 


Because of the prohibitive computational demands of three-dimensional 
(3-D) radiative calculations, most atmospheric dynamical simulations use one- 
dimensional (1-D) radiation models to calculate solar heating. However, several 
studies have indicated that the inherently multidimensional nature of cloud 
radiative processes can cause large errors in 1-D calculations (e.g., O'Hirok and 
Gauthier 1998, 2005; Di Giuseppe and Tompkins 2003, 2005; Pincus et al. 2005). 
Complementing the detailed case studies, 2-D radiative calculations using the 4 km- 
resolution clouds of a monthlong global simulation also showed some significant 
departures from 1-D theory (Cole et al. 2005). But even though such studies 
provided many valuable insights, it remains unclear how important the 
multidimensional nature of solar radiative processes is for atmospheric simulations. 

The US Department of Energy ARM Climate Research Facility (ACRF) offers 
excellent opportunities to address this issue by providing long-term detailed cloud 
observations at several sites. At these sites, ground-based vertically pointing 
instruments provide 2-D vertical cross-sections through the clouds drifting aloft. By 
comparing the results of 1-D and 2-D radiation simulations for clouds observed at a 
tropical, mid-latitude, and arctic location, this paper provides some initial estimates 
on the typical magnitude of 2-D radiative effects. Because the calculations don’t 
include cross-wind cloud variability, the presented 2-D effects can be considered as 
conservative, lower-bound estimates for the full 3-D effects (Pincus et al. 2005). 
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2. Dataset and methodology 


We analyze three years of cloud data at ACRF's NSA (Northern Slope of 
Alaska) and SGP (Southern Great Plains) sites, from 2005-2007 and 1999-2001, 
respectively. At the TWP (Tropical Western Pacific) site we analyze the two years 
(2003-2004) for which Microbase data is available for the full year. The ACRF 
Microbase product combines millimeter-wavelength cloud radar, micropulse lidar, 
and microwave radiometer data to estimate cloud liquid and ice water content and 
particle size at 10 s temporal and 45 m vertical resolution. These time-dependent 
vertical profiles are converted into 2-D spatial structures using each cloud layer’s 
mean wind speed from the ACRF Mergesonde product and the frozen turbulence 
assumption (e.g., Barker et al. 2004). Though the actual resolution of the obtained 2- 
D fields varies with the wind speed (median values: 86 m at NSA, 141 m at SGP, and 
74 m at TWP), for the radiative calculations we resample all data to a uniform 50 m 
resolution. 

While the ACRF Microbase and Mergesonde products are certainly not 
perfect and are influenced by both instrument uncertainties and retrieval algorithm 
limitations, they arguably represent the current state-of-the-art for comprehensive, 
long-term, and detailed cloud structure datasets. This analysis considers all daytime 
data, but excludes 10 km wide swaths around data gaps and around data flagged 
with a low quality control flag. 

We simulate radiative transfer through the observed clouds using a forward 
Monte Carlo model that was tested in the I3RC project (Cahalan et al., 2005) and 
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through comparisons to broadband SBDART (Ricchiazzi et al. 1998) calculations. 
The 2-D simulations assume no variability in the cross-wind direction and — to 
capture as much of the 3-D radiative process as possible — they assume that the 
solar azimuth is parallel to the wind (e.g., Varnai and Marshak 2003, Pincus et al., 
2005). We perform two 2-D simulations for each scene, with the sun on the 
downwind and upwind sides, respectively. 

Because the goal is to obtain robust statistics on the difference between 1-D 
and 2-D results — as opposed to obtaining accurate results for individual clouds — 
the simulations can use both relatively few photons (up to about 3000 for each 50 m 
column), and moderately complex radiative characteristics. Minor imperfections are 
not expected to influence the obtained statistics because the 1-D and 2-D 
simulations use the same Monte Carlo radiative transfer model and setup. 

The Monte Carlo model calculates gaseous absorption in the 0. 2-5.0 pm 
range using the correlated-k method, with coefficients from SBDART (Ricchiazzi et 
al. 1998). The SBDART tropical atmospheric profile is used for the TWP site, and the 
mid-latitude summer (winter) and sub-arctic summer (winter) profiles are used for 
the SGP and NSA sites during the April-September (October-March) periods, 
respectively. Cloud particle scattering and absorption parameters are from Mie 
calculations for liquid droplets, and from the publicly available database of Baum 
(http://www.ssec.wisc.edu/~baum/Cirrus/Solar_Spectral_Models.html) for ice 
crystals. Rayleigh scattering is considered, but aerosol effects are not included. 

Since Barker and Davies (1992) shows that the angular pattern of surface 
reflection has only a modest influence on 3-D radiative effects, all simulations use 
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Lambertian surfaces with 5 nm-resolution spectral albedos. The sea surface albedos 
used at the TWP site are from the database described in Jin et al. (2004]. The used 
albedo values depend on the local mean cloud optical thickness and solar elevation, 
and assume a 5 m/s wind speed. Surface albedos for the SGP and NSA sites are 
specified by the SBDART model's "vegetation", "sand", and "snow" surface types. At 
SGP, the surface type is selected daily based on human observer reports about 
surface conditions (documented in ACRF’s Surflog product], with the "sand" type 
used for days without snow or green vegetation. At NSA, where no such data is 
available, "vegetation" is used from June to September, and "snow" for the rest of the 
year. These dates were selected based on snow cover data from the Rutgers 
University Global Snow Lab ( http : / /climate.rutgers.edu /snowcover /] . 

3. Results 

We examine the influence of 2-D radiative effects by comparing the results of 
1-D and 2-D radiative simulations. Table 1 presents the influence of 2-D effects on 
multiyear 24-hour average solar radiative fluxes (including cloud-free periods and 
nighttime as well]. The table shows that 2-D effects are stronger at lower latitudes, 
due to factors such as 

• more intense solar illumination, 

• stronger convection creating more heterogeneous clouds, 
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• more effective channeling process for high sun, guiding photons from opaque to 
thin regions where they can reach the surface easier (e.g., Davis and Marshak 
2001 ). 

We note that the 2-D effects in Table 1 may be considered as conservative, 
lower-bound estimates for the full 3-D radiative effects, because earlier case studies 
found 3-D effects to be about 30% stronger than 2-D effects (e.g., Pincus et al. 2005). 
The difference arises because 2-D calculations don't consider the cross-wind cloud 
variability (perpendicular to the solar azimuth) that further reduces cloud reflection 
and enhances surface absorption through additional channeling. We also note that 
in coarse-resolution cloud simulations, the errors of 1-D albedo and surface 
absorption calculations are further increased by the plane-parallel bias caused by 
unresolved small-scale cloud variability (e.g., Cahalan et al. 1994). As a final note, we 
also mention that for reflective surfaces, the influence of 2-D effects is somewhat 
larger on downwelling fluxes at the surface than on surface absorption. 

Figure 1 shows the way average 2-D effects depend on solar zenith angle. 2-D 
effects being largest for high sun means that they are strongest near noon and 
during the summer. The figure also shows that for relatively high sun, 2-D effects 
are nearly as strong at NSA as at SGP, which implies that Alaska clouds in the 
summer and near-noon are almost as heterogeneous as Oklahoma clouds in the 
spring or fall, or in the morning or afternoon. 

The figure also shows that 2-D effects increase atmospheric absorption even 
for high sun. This may be explained by cloud side illumination increasing the 
sunlight intercepted by clouds, and also by the trapping of sunlight in thicker cloud 
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portions where it is more likely to get absorbed (e.g., O'Hirok and Gauthier 1998; 
Varnai and Davies 1999). 

A somewhat unexpected feature of Figure 1 is a behavior opposite to the 
dominant 2-D effect discussed in most earlier case studies (e.g., McKee and Cox, 
1974) where, for low sun, cloud sides intercepting extra sunlight enhanced 
reflectivity and reduced transmission to the surface. Figure 1 indicates that 2-D 
effects often enhance surface absorption even for oblique sun. As Figure 2 
illustrates, this occurs when the incoming direct sunlight slips obliquely under an 
extensive cloud and, after reflection from the surface, gets trapped by the cloud 
above (O’Hirok and Gauthier 1998). 

Naturally, 2-D effects also depend on cloud type. We examine this issue using 
the ACRF cloud classification product (Cldclass) available at the SGP site (Wang and 
Sassen, 2001). Figure 3 shows that 2-D effects tend to be strongest for the more 
heterogeneous, convective cloud types. The main exception is surface absorption 
near deep convective clouds, because these highly reflective clouds greatly reduce 
surface absorption in both 1-D and 2-D simulations. 

While average values can be informative, detailed histograms can also help in 
evaluating the importance of 2-D effects. Figure 4 shows that local 2-D effects are 
often much larger than the average values discussed above. For example at TWP, 
2-D effects change the total cloud absorption of 1 km-size columns by more than 
20% in about a third of cases, and they change the total surface absorption of 
10 km-size areas by more than 50 W/m 2 in about a quarter of cases. 
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4. Summary 


This study presents multiyear statistics on the influence of 2-D shortwave 
radiative effects that are not considered in the 1-D radiation calculations used by 
most atmospheric dynamical models. The influence of 2-D effects on solar radiative 
energy budget calculations is presented for a tropical, mid-latitude, and arctic site of 
the US Department of Energy ARM Climate Research Facility. The results show that: 

• 2-D effects increase 24-hour average total solar absorption by about 4.1 W/m 2 at 
the TWP site in Papua-New Guinea, by 1.2 W/m 2 at the SGP site in Oklahoma, 
and by 0.3 W/m 2 at the NSA site in Alaska. 

• 2-D effects are often much larger than these average values, especially for high 
sun and for convective clouds. Such variations in 2-D effects can change the 
spatial and temporal distribution of solar radiative heating. 

• 2-D effects often increase solar heating even for oblique sun. This occurs when 
the incoming direct sunlight slips obliquely under an extensive cloud, and the 
cloud subsequently intercepts much of the light reflected from the surface. 

The presented 2-D effects can be considered as conservative, lower-bound 
estimates for the influence of horizontal photon transport not included in 1-D 
radiation calculations, because 2-D calculations don’t consider cross-wind cloud 
variability. In earlier case studies the full 3-D effects exceeded 2-D effects by about 
30% (e.g., Pincus et al. 2005). We note that the upcoming installation of scanning 
radars at ACRF sites will offer new opportunities for fully 3-D studies. 
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The results imply that considering horizontal photon transport can greatly 
improve both the interpretation of shortwave radiative flux measurements and the 
solar heating calculations in atmospheric dynamical models. This underscores the 
need for fast radiation calculation methods that can allow dynamical models to 
consider the inherently three-dimensional nature of shortwave radiative processes. 
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Tables 


Table 1. Multiyear 24-hour average 2-D effects on shortwave radiative fluxes, 
including nighttime and clear areas as well. The values indicate the difference 
between 2-D and 1-D fluxes simulated at full resolution. Uncertainties are standard 
errors based on the spread of results when the full dataset is divided into 25 
subsets. These standard errors likely underestimate actual uncertainties, because 
they assume that the 25 subsets are independent from each other. 


Site 

TOA reflected flux 
CW/m 2 ) 

Atmospheric 
absorption (W/m 2 ) 

Surface absorption 
CW/m 2 ) 

TWP 

-4.10 ±0.21 

1.47 ±0.03 

2.63 ± 0.17 

SGP 

-1.21 ±0.05 

0.52 ± 0.02 

0.68 ± 0.05 

NSA 

-0.28 ± 0.04 

0.25 ± 0.01 

0.02± 0.04 


13 




Figures 



Flux absorbed at the surface 



Figure 1. Average 2-D effects on radiative fluxes, as a function of solar zenith angle, 
(a) Reflected flux; (b) Flux absorbed in the atmosphere; (c) Flux absorbed at the 
surface. Error bars in this and subsequent figures show standard errors based on 
the spread of results when the full dataset is divided into 25 subsets. These standard 
errors likely underestimate actual uncertainties, because they assume that the 25 
subsets are independent from each other. 
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Distance (km) 


Figure 2. Illustration of the photon trapping process that enhances surface 
absorption for oblique sun. (top] Illustration using an observed field of volume 
extinction coefficients; (bottom] Influence of 2-D effects on simulated surface 
absorption values. 2-D effects enhance the scene average surface absorption by 
9.6 W/m 2 . The solar zenith angle is 69°. 
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Figure 3. Average 2-D radiative effects for various cloud types at the SGP site, 
(a) Relative change in cloud absorption; (b) Absolute change in daytime average 
surface absorption within 10 km from clouds of each type (W/m 2 ). 
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Figure 4. Cumulative histogram of 2-D effects on radiative fluxes. Each curve is for a 


different year, (a) Relative change in the total absorption of 1 km-size columns; 


(b) Absolute change in the daytime average surface absorption of all 10 km-size 


areas, including clear regions. 
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